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Abstract. We consider the process 77 H1H2 where Hi and H2 are either mesons or baryons. The 
experimental findings for such quantities as the pp and KsKs differential cross sections, in the energy 
range currently probed, are found often to be in disparity with the scaling behaviour expected from hard 
constituent scattering. We discuss the long-distance pole-resonance contribution in understanding the 
origin of these phenomena, as well as the amplitude relations governing the short-distance contribution 
which we model as a scaling contribution. When considering the latter, we argue that the difference 
found for the KsKs and the K'^ K~ integrated cross sections can be attributed to the s-channel isovector 
component. This corresponds to the pui —> a subprocess in the VMD (vector-meson-dominance) language. 
The ratio of the two cross sections is enhanced by the suppression of the component, and is hence 
constrained. We give similar constraints to a number of other hadron pair production channels. After 
writing down the scaling and pole-resonance contributions accordingly, the direct summation of the two 
contributions is found to reproduce some salient features of the pp and K~ data. 

PACS. ll.30.Hv Flavour symmetries - 12.40.-y Other models for strong interactions - 12.40.Nn Regge 
theory, duality, absorptive/optical models - 12.40.Vv Vector-meson dominance - 13. 66. Be Hadron pro- 
duction in e~e^ interactions 



1 Introduction 

We consider the exclusive pair production process 77 — > 
H1H2 where Hi and H2 are either mesons or baryons. We 
consider the energy region not too far from the threshold, 
for example the centre-of-mass energy W^^ < 4 GeV. 

It has been found that some of the results of recent 
large-statistics measurement of these processes, for exam- 
ple at Belle jT|[2] , are difficult to reconcile with theoretical 
thinking based on hard constituent scattering. For exam- 
ple, there is a violation of the scaling behaviour expected 
from the naive quark-counting rule 

da cos 9* , , 



where K is the number of 'elementary' fields taking part 
in the interaction. For instance, _fi' = l + l + 2 + 2 = 6 
for 77 — > TTTT. s and t are the usual Mandelstam variables. 
After integration over a constant cos Q* interval where Q* 
is the polar angle of scattering in the centre-of-mass frame, 
this yields: 

a cx Xj (mesons), (baryons). (2) 
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This scaling is expected to hold when the constituent- 
level hard subprocess approximately has a tree-level per- 
turbative description. As seen in fig. [l] this fails for 77 
pp mmSlE] in the measured energy range. It seems to 
work for 77 ^ K+K- pIlTllS] above W^^ « 2.4 GeV but, 
disturbingly, seems to fail for 77 — > KsKs 9l in the mea- 
sured energy range of up to W^^ « 4 GeV. The latter cross 
section drops faster with W^^^, and this finding, together 
with the large ratio between the two cross sections, is dif- 
ficult to explain in frameworks based on hard constituent 
scattering. On the other hand, in the large-energy limit 
of the measured range, it has been found that the calcu- 
lations of ref. [ini, based on leading-term QCD and wave 
functions following from the QCD sum rules, can accom- 
modate the value found experimentally. 

Our main region of interest in this paper is below this 
energy range, where pcrturbative description is insuSicient 
to account for the prominent features of the data. We 
are interested in the participation of alternative dynamics, 
that are more long-distance in nature, and are more ap- 
propriate to describing the observed distributions. At the 
same time, we are also interested in modelling the short- 
distance contribution with constraints from the amplitude 
factorization considerations. 

As a starting-point, let us consider VMD (vector- meson 
dominance). Here, the photon is interpreted as a quark- 
antiquark object, so that the exponent K in eqn. ([!]) is 
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give the physical / and /' mesons. We are not necessarily 
adopting the s-channel resonance picture, and / and a are, 
for now, merely a label of the s-channel flavour structure. 



OAAA 




(a) 



(b) 



(c) 



Fig. 2. The quark- line diagrams for meson-pair production. 
Similar diagrams can be drawn for the baryon-pair case. Dia- 
grams (a) and (b) have s-channel representation whereas the 
4-quark mode (c) does not. 



Fig. 1. 77 — > pp cross section [T] versus centre-of-mass energy 
at VENUS g], CLEO and Belle [1] in the central region, 
defined by | cos8*\ < 0.6. The vertical error-bars on the Belle 
data are due to the statistical error in the event and the Monte 
Carlo samples only. Experimental data is compared against 
three theoretical calculations, as described in ref. [T|. Figure 
reproduced with kind permission of the authors of ref. [l] . 

modified. Corresponding to eqn. ([2]), we would have: 

cr cx (mesons), l/s"^ (baryons). (3) 

Although the applicability of the quark-counting rule to 
the VMD picture should not be taken for granted, this 
indicates that the fall in the cross section with the centre- 
of-mass energy would be more rapid than is expected from 
eqn. ([1]). It is a curious finding that for central events, 
defined by | cos 6**1 < 0.6, the cross sections measured at 
Belle for 77 — > KsKs 9^ and pp [I] go as ~ W^-/ and 

^ W^^^'^'^^^^ respectively for some regions of W^^ away 
from the resonance region. As the exponent is sensitive to 
the cut on | cos 6** | , it is possible that this agreement with 
eqn. ([3]) is accidental. We note nevertheless that a result 
of the form above can be derived from the consideration 
of the Sudakov form-factor effects [TT[|T2| . 

Encouraged by this finding, we go on to consider the 
factorization of the scattering amplitude into the produc- 
tion and decay parts. The production subprocess is dom- 
inated, in case of ideal mixing between uj and 0, by: 

p°p\LOLO^Ud. (4) 

pV^P°^a, (5) 
H^fs. (6) 

In the above, fud stands for the {uu + dd)/^/2 state. After 
relaxing the condition of ideal mixing, fud and fs nrix to 



The quark- line diagrams are shown in fig. El The production- 
decay factorization assumption implicitly corresponds to 
the choice of 'handbag' diagrams (a) and (b) of fig. O ne- 
glecting the 4-quark intermediate state of 'cat's-ears' dia- 
gram (c). When considering long-distance dynamics, due 
to the degeneracy of the meson trajectories, we shall argue 
later that this approximation is acceptable for meson-pair 
production but not for baryon-pair production. On the 
other hand, while calculating the short-distance contribu- 
tion, the relative size of the cat's-ears diagram is sensitive 
to the factorization scheme adopted. Here we limit our- 
selves to remarking that in the findings of ref. |13| , and 
for the energy range currently probed by experiment, their 
so-named 'handbag' approach is appropriate for describing 
some features of the data at hand. For instance, the angu- 
lar distribution for the KsKs process, as seen in ref. 
is in better agreement than the more traditional approach 
of ref. [10 , although the overall normalization is not well 
understood in either of the two approaches. 

We have also neglected the OZI-suppressed channels 
such as — > /. The decay part can be expressed simi- 
larly. 

We see immediately that any difference between the 
K^K~ and K'^K'^ = K^Kl + KsKs cross sections must 
be due to the simultaneous presence of the isoscalar ///' 
and isovector a = {uu — dd) / components. It is also not 
difficult to see that the large ratio between K'^K~ and 
KsKs processes can be increased by the suppression of 

(j^ff) Is- 

We proceed by modelling the short-distance piece as a 
scaling contribution obeying eqns. (|ll2p . The long-distance 
piece includes the resonances of eqns. (|4l)-([6]). These are 
related to the i-channel pole picture by duality, so that 
we can also model them as Regge amplitudes [l4l[T5] . The 
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signature term in the latter case would then represent the 
'cat's-ears' contributions. 

We find that the direct summation of the two ampli- 
tudes reproduces some of the salient features of the K^K^ 
and pp data, although the description of the angular dis- 
tribution in the former case is poor. 

The modelling of the full amplitude as the sum of 
the two contributions at first sight may seem to suffer 
from the problem of double counting. However, we find 
numerically that adopting the alternative approach of a 
form factor that interpolates between the two amplitude 
in general cannot yield these results. An intuitive expla- 
nation would be that the finite-time short-distance effects 
and infinite-time Regge-pole dynamics have little over- 
lap. More explicitly, short-distance amplitudes involve a 
finite and small number of intermediate particles that in- 
dividually carry large virtuality, and are independent from 
the long-distance amplitudes where virtuality is assigned, 
in the parton language, collectively to the intermediate 
state partons. Hence in our understanding, any discrep- 
ancy with the data that arise are due to our ignorance of 
low-energy dynamics, rather than being due to some form 
of double-counting in this summation procedure. In any 
case, such approach is not new. For a well-known exam- 
ple involving the summation of short-distance and long- 
distance amplitudes, see ref. [IB]. 

This paper is organized as follows. In sec.[2l we present 
the relevant Clebsch-Gordan coefficients for kaons and for 
other mesons and baryons. In sec. [31 we consider the effect 
of pole (-resonance) dynamics. The conclusions are stated 
at the end. 



2 SU{3) analysis 

We decompose the (2 2) amplitudes into the s-channel 
production part and the decay part, so that, for example, 
the s-channel scalar part of the amplitudes is given by: 

= IV1IV29V1V2S X gsHiH2 X F{S). (7) 

g are proportional to the SU{3) flavour Clebsch-Gordan 
coefficients, whereas the dynamics is contained in the func- 
tion F{S). 7y are the photon-F coupfing constants satis- 
fying: 



7p0 



7t 



3 : 1 



-V2. 



(8) 



We shall adopt the ratio 3 : 1 for the p° and uj couplings 
later on, but modify the (f> coupling by a suppression factor 
VS- It should be understood that this suppression factor 
is introduced in order to suppress /' contribution. The 
photon coupling to 4> is not changed. 

We note that the pp data [T] shows clear indication 
of the presence of a pseudoscalar rjc resonance, and the 
same peak is also present in the KK data The pseu- 
doscalar contribution can be included without modifying 
the structure of the formalism. 



2.1 Production subprocess 

u, d and s quarks are organized into a flavour triplet struc- 
ture as: 




3* = = (u J 



(9) 



Mesons are in 3 3* and baryons are in 3 (8) 3 (8) 3. 

We define the nonet 1~, 0~, and 0+ mesons by V^, , 
and S*^, respectively. These are constructed explicitly as: 



uu ud us 
3* = I du dd ds 
su sd ss 



p +0; cos 0v — sin ffy + 

— —p^-\-uj COS Oy —<f> sin 6v 

K*- K*o 



(10) 



K*+ 

K*o ] (11) 

Lu sm fJv + (p cos Ov , 



V6 



Vs 



V3 



K+ 



TT 

K 



V6 



V3' 



na 
Oi, 



UU Kq ' 

Gq dd Kq 



(13) 

mixing angle Oy in the above. 



We have defined the lu — 
For ideal mixing, 9v = 0- 

The Clebsch-Gordan coefficients for the production sub- 
processes of eqns. (H])-® are calculated by the contribu- 
tion of the diagonal, a — b — c, part of the quantity: 



(14) 



This particular notation implies spin-0 state in the s- 
channel, but the structure of the expression is general to 
any even-spin positive parity states. The / — /' mixing 
angle 9s is defined by: 

flo = (uu — dd)/V2, (15) 
fo ^ cosOsiuu + dd)/V2 + sinOsiss), (16) 
/o = -sin9s{uu + dd)/V2 + cosdsiss). (17) 

For later use, we define 9p analogously, which describes 
T] — rj' mixing. 

The Clebsch-Gordan coefficients can now be calcu- 
lated, and these are listed in tab. [TJ We show the general 
case as well as the two cases of ideal mixing, correspond- 
ing to 0y = and to Oy = Os = 0. In reality, 6*1/ = is a 
good approximation, although ^5 = is doubtful for the 
low-lying resonances. The terms that vanish in the ideal 
mixing case are OZI-suppressed. 

One consequence of the OZI suppression is that when 
we consider the full production process, namely 77 — > 
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vertex 


general mixing 


6v ^0 


Sv- = = 


P"p" - fo 


cos 


cos 6ls/V2 




P°p' - /o 


-sin 6*5/72 


- sin 6s /\/2 





p°LO,U!p'^ ao 


cos 9v /\f2 


l/%/2 


l/%/2 


p°(j), cf)p° ~* ao 


-sin ev/\/2 








UlUJ — > fo 


cos^ 9v cos 6s 1 \/2 + sin"^ Qv sin Qs 


cos 6s /V2 


I/V2 


LOUJ ^ fo 


— cos^ Qv sin Qs / \/2 + sin^ Oy cos 6s 


- sin 6s /V2 





LU4> /() 


sin dv cos dv{— cos 6s + sin 








uj4>-* fo 


sin 6v cos (sin 63/ \/2 + cos 6's ) 








cj)(j) fo 


sin"^ 6v cos 0s /\/2 + cos'^ 6v sin 


sin 6 s 





H fo 


— sin^ 6v sin 6*5 /•/2 + cos^ 6*^ cos 6s 


cos 


1 



Table 1. The VVS coupling coefficients. The overall symmetry factor 2 has been suppressed. The vertexes not listed here are 
forbidden by isospin conservation, so that we have: gpOpOag = 9p0^fa = ffpO^,/^ = 5pO.^>/o = 5pO0/^ = P^^^o = S'L.^ao = 50<#.ao = 0. 
We list the general case as well as the 'ideal mixing' case corresponding to dy ~ and to 6v ~ 6s = 0. 



V1V2 with the vector mesons in the final state, the channels 
which necessarily involve OZI-suppressed interaction are 
suppressed, so that we expect: 



f (77 ~* P^p'^)j '^{11 ~* 

> (7(77 — > p'^Ijj) > (T{"f"f ■ 

» 17(77 ~^ p^4 



,cr(77 ^ w0).(18) 



The reasoning goes as follows. Let us first emphasize that 
this is for the entire s-channel production process as given 
by the factorization of eqn. The photons first couple 
to the appropriate vector boson, which fuse together into a 
(scalar) resonance, then finally decay into the states given 
above. In the ideal mixing case, p°p° and both come 
from / decay with equal strength, so that and ujuj are 
approximately equal, p^uj can come from a, but a produc- 
tion is slightly smaller than / because of the ratio of p° 
and oj contents of the photon, (jxj) is suppressed because /' 
is less abundant than /, that is, p^ f is the dominant 
resonance production subprocess. 



2.2 Meson pair production 

We now turn to the decay subprocess. Let us first consider 
the production of pseudo-scalar mesons. These are in the 
nonet representation of eqn. IT^ . 

The relevant coefficients can be obtained by the diag- 
onal, a — h, part of: 



PcPb)- 



(19) 



We list them in tab. H 

Again, those modes that are OZI suppressed are ac- 
companied by factor smOs or smOp. However, since Bp is 
now considerably large, ~ —39 degrees [T7], the suppres- 
sion factor is only moderate. 9s is also considerably large 
for the spin-0 bosons T8 although possibly not for the 
higher-spin excitations of a, / and /'. 

We see that, after including the identical particle factor 
of 1/2 for the tt^tt" cross section: 



cr(77 ■ 



-) = 2cr(77 -> tt'V") 



(20) 



regardless of the mixing angles or production dynamics. 
We also see that the channels: 



77 ■ 



(21) 



can only proceed via s-channel isovector a, so that the ob- 
servation of these processes would be interesting to con- 
firm the presence of the isovector channel. The magnitude 
of these channels is related to the difference of the kaon 
amplitudes, i.e.: 



^(77 K+K-) - A{-f-f ^ K"K0) 



^(77 — > TT 77) + LA(77 — > 



(22) 



In particular, when the K^K 



cross section dominates 
over K^K'^, and when W^^ is sufficiently above the tt^ti' 



threshold, we obtain: 
(7(77 K^K^) K, (7(77 



7r°77) -|-(7(77 -> 7r°77')- (23) 



This should be tested experimentally. However, we shall 
show later that the ratio between the K~^K~ and K'^K'^ 
amplitudes is expected to be at most 4, and hence omitting 
the K'^K^ amplitude contribution is not a good approx- 
imation. One possible assumption would be that the two 
cross sections only differ by a real constant factor. The 
cross section would then be scaled by: 

(7(77 ^ K+K-) — > 



a{j-f K+K-) X 1 



(7(77 ■ 



1(24) 



a{-fj ^ K+K-) 



Similarly, for the sum of the K^K 
tudes, we find that: 



^(77 K+K ) + ^(77 
~ -4(77 — > r/77) + ^(77 - 



and K^K'^ ampli- 



V V ) 



(25) 



However, this equality would be difficult to test experi- 
mentally, unless it is found, for instance, that rji] domi- 
nates over rj'rj' when sufficiently above the threshold. In 
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final state oq /o 7o~ 








COS 9sV^ 




K+R- 


1 


^+ sines 


cos 0s 


TV TT 





cos 6sV2 


-sin6isV2 


K°7{0 


1 

V2 


^+ sines 


+ coses 



Tv°rj, r?7r" cos 6lp\/2 

TT°ri',ri'-K° -sin6'p^/2 



777? 2 (^cos^ e/p i^2|s + giji^ 6lp sin ^s j 2 cos^ 6lp + sin^ 6»p cos ^s 

7???' 2smepcosep{~-^2^ + sines) 2 sin Sp cos 6lp(2i^ + cos 6's) 

T?'?7' 2 f sin^ 6)p + cos^ Sp sin 615) 2 f - sin^ Sp + cos^ Sp cos 63 



Table 2. The S'PP coupling coefficients. 



this case, the relation would reduce to: 



1 



(7(77 K+K ) w -cr(77 -f rjT]) 



(26) 



The i factor on the right-hand side comes from the com- 
bination of the factor 2 in eqn. (j25|) and the factor ^ for 
identical particle production. Again, we can adjust for the 
error in neglecting the sub-leading amplitudes by a scal- 
ing similar to eqn. (|24p but with the minus sign in the 
brackets replaced by a plus sign. 

The above argument holds regardless of the production 
subprocess VV S. Let us now consider the inclusion of 
the production subprocess. This is obtained by referring to 
eqn. ([7]). We adopt ideal mixing for the vector mesons so 
that 9v = 0. The result, for 77 K+ R- , K^W ,tt+tt- , 
and 7r°7r° are: 

y^(77 ^ K+K~) 

/ cos6's ■ a \ ^°^^S , S ■ Q \ 

^Ssin^s S \ 

;= cos Us 



sin 6s n \ I ^ '^i" " 
— cos Os ] [ - — 1= 



(27) 



/ cos6'5 ■ a \ ^°^^S , S ■ Q \ 



^(77 TT+vr^) 

- F(/o) cos 0s (^^ cos 0s + ^ V2 sin 

+F{f;,) sines (^^sin^s- ^V2cos0s) , 
^(77 ^ ttV) 

= F(/o)cos0s (^^cos0s + ^V2sin0s) 



"Ssin^s S 
^ ;= cos 



SUs 



(29) 



+F(/^)sin0s (^^sin^s- ^V2cos0s) 



(30) 



5 is defined below eqn. ([5]). The above expressions simplify 
in the case 9s — 0. As mentioned earlier, this becomes 
acceptable for the higher-spin excitations of /, a, and /'. 
We then have: 

yl(77 -> K+K-) 

= ^[5F{fo) + 2SFif;,)+3F{ao)], (31) 
^(77 ^ i^°!^) 

= ^ [5i^(/o) + 2<5F(/^) - 3F{ao)] , (32) 



^(77->vr+7r-) = -i^(/o), 
^(77 -> ^°^°) = ^^(/o). 



(33) 
(34) 



We may simplify further by the approximation -F(/o) = 
F{ao) and absorbing the difference between -F'(/o) and 
F(/o) in the coefficient 6. Under these conditions, i.e., 
F{a) — F{f) — F{f') for arbitrary 5, it turns out that 
all dependence on 9s cancels automatically, so that ideal 
mixing becomes a redundant assumption. We obtain, at 
the amplitude level: 

K+R- : is:"ifO : tt+tt" : 7r°7r° = A + 5 : 1 + 5 : 5 : 5. (35) 



In particular, for the ratio of the 77 
cross sections, we have: 



K+R- and RsRs 



(28) cr(77 _> R+R-) _ 2^^(77 K+R-) 
(7(77 RsRs) ~ cr(77 ^ i^o^o) 



1 



(36) 



The conventional charge-counting argument [13j corresponds 
to (5 = 1 and therefore gives 25/2. On the other hand, the 
complete suppression of /' gives rise to 32. Since physi- 
cally we expect 5 to be positive definite, we obtain the 
following inequality: 



12.5 < 



cr(77 -> R+R-) 
(7(77 RsKs) 



< 32. 



(37) 
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The upper limit of this equation seems to be satisfied by 
the currently available data [9], within the statistical er- 
rors. The closeness of the observed ratio at high energy 
with the limiting value 32 indicates the suppression of the 
strangeness coupling in this region. The lower limit seems 
to be violated at lower energies [9] . One possible reason is 
that, in this region, there is dominant contribution from 
/' corresponding to the possibility 5 > 1, but another pos- 
sible reason is that, as we shall argue later, long-distance 
interaction tends to respect isospin invariance, and so the 
a coupling is suppressed. 

Another quantity that may be interesting is the ratio 
of the K~^K~ and tt+tt" cross sections. We obtain: 



cr(77 ^ K+K~) 
(7(77 — > 7r+7r~) 



(38) 



We therefore expect that sufficiently above the threshold, 
the K^K'^ cross section is smaller than the tt+tt" cross 
section by up to a factor of 16/25. 

If the suppression factor S is indicative only of the res- 
onance structure denoted by F{f'), and not the couphng, 
of /', it then follows that the 5 obtained by the compar- 
ison of these cross sections should be approximately the 
same as the S obtained by the comparison of the K~ 
and KsKs in the same kinematic range, i.e., the same 
W-y^ and cos9* intervals. This should be tested experi- 
mentally. 



2.3 Baryon pair production 

Baryons belong to the 3 (g) 3 (g) 3 representation, which is 
decomposed as: 

3(8) 3(g) 3 = 10s eSMS © 8ma © 1a- (39) 

Out of these, phenomenologically the most relevant are the 
1/2+ baryons in the octet representation. The subscripts 
MS and MA stand for mixed symmetric and mixed an- 
tisymmetric, respectively. We write them as: 



8ms — 5(a,6)ci 
8mA — B[a.b]c- 



(40) 
(41) 



Curly brackets in the subscript represent the symmetric 
sum and square brackets represent the antisymmetric sum, 
so that: 

B(a,b)c = B(^i„a)c, B[aM]c = —B[i,^a\c- (42) 

Furthermore, -B[a,6]c ^Iso satisfy the Jacobi identity: 

B[a,b]c + S[(,^c]a + B[c,a]b = 0. (43) 

The symmetric and antisymmetric representations contain 
the same physical states, and they are related by: 



1 



-Of) = e ^{d.b)c — -D[c,d]67 



(44) 



where e'"^'' is the Levi-Civita tensor with the convention 
£123 = i_ is the octet m; 
content of B[a.b]c explicitly: 



gi23 _ -|^^ j^a |.]-^g octet matrix. We first write down the 



B[ia]i = P, -S[i^2]2 = n, -B[i,2]3 = -"^A/VG, 



B,2^3U^S^/V2 + A/V6.B, 



The octet matrix B^ is then: 



2,3]2 



[2,3]3 



-H(45) 



B^ 



'uyV2 + A/V6 S+ p 

S- -i;0/\/2 + yl/%/6 n 



(46) 

There are three possibilities for evaluating the baryon- 
baryon-meson SBB coupling. We can work in terms of 
B(a,b)ci ^[a,6]c Or B^ . Hcrc we adopt the following nota- 
tion: 

iaB''''l'^i?[,,,]fcA4^ + /?B'"''l'S[b,,]rfM^ (47) 



a and /3 are adjustable parameters satisfying [19:: 



a ~ 513. 



(48) 



This relation comes from the approximate flavour SU{3) 
symmetry for baryon- meson strong coupling constants. 
On the other hand, for the naive charge-counting argu- 
ment to work, we need to impose a = /?. This point will 
be demonstrated by an example later. 

We can now tabulate the relevant coupling constants 
in terms of a and /3, and these are listed in tab. [3l 

We can consider forming equalities similar to eqns. ([22|) 
and ([^5]) . In particular, we can make use of the ratios of 
flQ couplings: 

.4(77 ~+ pp) ~ .4(77 nn) : 



^(77 S+S+) - 
^(77 -> H"^) - 
^(77 ^ S"A) 

^ a : a + f3 : P 



a- (3 



2V3 



(49) 



However, these will be difficult to verify experimentally. 
This is partly because some final states, for example nn, 
are difficult to measure, and partly because we do not 
expect in any of the pairs of reactions above that either 
of the two amplitudes would become sufficiently dominant 
over the other that the other can be neglected. 

It is hence more helpful to make an estimation analo- 
gous to eqn. psp . For example, the pp amplitude is given, 
for Ov = 0, by: 

^(77 pp) 

= ^(a + 2/3) [cos^ esF{So) + Jsin^ OsFUl,)] 



:F{ao) 



(50) 
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final state oo Jo 

pp a/V2 (a + 2/3) cosOs/V^ -(a + 2/3) sin 63 /V^ 

n n -a/V2 (a + 2f3) cosOs /V^ -(a + 2/3) sin6ls/\/2 

IJ+E+ ^ (a + l3)cos9s/V^ + Psines ~{a + f3) sinOs /y^ + f3 cosOs 

(q + /3) cos 6is/ V2 + /3 sin 6is -(a + /3) sin 6*5/^2 + /3 cos 6is 

IJ-E^ (Q + /3)cos6ls/%/2 + /3sin6is -(a + /3) sine»s/x/2 + /3cos6is 

AA ^^tMcos^s + ^^sin^s -2±^ sinOs + ^ cosOs 

S°HO /3/%/2 (q + /3) sin e»s + /3 cos 6*5/^2 (a + /3) cos6is - /3sine»s/V2 

-13 /V2 {a + I3)sm6s + Pcos6s/V2 (a + /3) cos 0s - /3 sin 6*5/72 

i:°A,ASO ^ 

Table 3. The 5-6-6 coupling coefficients. 



As before, by taking ds = 0, F{fo) = -F(ao) and absorbing 
the difference between F{fo) and -F(/o) into the coefhcient 
S, we arrive at: 



A{j-f-*pp) « -(4a + 5^) 



(51) 



We repeat the same exercise for the other production 
modes and obtain results listed in tab. [H 

The statement made earlier about naive charge-counting 
in the case a = (3 can now be demonstrated explicitly. For 
example, the ratio of the pp and nn amplitudes is given 
by: 

8 + 10 3 12 + 22 + 2^ 

^- (52) 
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1^ 



1^ 



The ratios of cross sections sufficiently above the thresh- 
old region are given by the square of the coefficients listed 
in tab.m so that, for instance: 



ia + 10/3 



o'(77 pp) ^ 
cr(77 ^ Z-OZ^) V5a + (5 + 2S)(3 



(53) 



This quantity comes out to be between ^ 2.4 and ~ 2.7 for 
a = 5 /3 an d < S < 1. The measurement of the S^E'^, 
or cross section would be particularly interest- 

ing because of the sensitivity to 6. From the discussion of 
pseudo-scalar meson pair production, we expect that 6 is 
small, so that this cross section would be suppressed by 
factor ^ (50/8)2 ^ compared to the pp cross section. 



3 Long- and short-distance dynamics 

Let us model the short-distance amplitude as a scaling 
contribution: 



A oc 



„4-K/2 



{t-Ml){u-Miy 



(54) 



K is as appearing in eqn. ([T|). We take the two masses 
Ml and M2 to be the corresponding hadron masses as 
opposed to, for instance, some appropriate quark masses. 
The contribution of these mass terms is negligible in any 
case in the region of interest. This expression gives rise to 



the angular distribution k, (1 — cos^ (?2)-2 -^^j^j^jj jg charac- 
teristic of single quark exchange, or more generally light 
particle exchange, in the ^-channel. This angular distri- 
bution would be valid at high-energy, and the scaling be- 
haviour of eqn. ([1]) implies that the angular distribution 
must remain the same at low-energy. 

To this, we add a long-distance pole(-resonance) con- 
tribution which has the Regge limiting form: 

At^r{l~ a{t)) (1 + r exp(-i7ra(t))) (a's)"^*) 
+ {t^u), (55) 

The linear trajectory is parametrized a{t) = q;(0) + a't as 
usual, with a' k, 0.9 GeV"^. ^ is the lowest spin of the 
trajectory, and r = ±1 is the signature. For baryons, the 
signature term is modified to |23| : 



1 + rexp {-iTT{a{t) - 1/2)) . 



(56) 



It is found [21] that a simple Regge expression similar 
to eqn. (|55|) yields characteristic behaviour in the central 
region cos 9* ~ that is in accordance with the 77 — > pp 
data just above the threshold. This may seem surprising at 
first sight, but is reasonable considering that the (Regge) 
pole amplitude and the resonance amplitude are related 
by (semi- local) duality i20j . After integrating over the res- 
onances, the behaviour of the two amplitudes is similar. In 
particular, this method works well in the case of pp since 
the cross section is a smooth function of W^^ and no trace 
of resonances is seen. For K^K~, the resonance struc- 
ture is still seen, so that we may, for instance, replace the 
Regge amplitude with a resonance-pole dual amplitude of 
the Veneziano model [221 • 



3.1 Baryon pair production 

Let us first consider pp production. In eqn. ([54]) . we set 
-fiT = 8. In eqn. ([55)) with the modification of eqn. ([5S)) . 



we set ^ = i. As for the trajectories, ref. [23], the lead- 
ing 5 = contributions, the N/A, have the following 
parametrization: 

(57) 
(58) 
(59) 
(60) 



iVa : 


a{t) 


= -0.34 + 0.99i, 


: 


ait) 


= -0.63 + 0.89i, 


Np : 


a{t) 


= + 0.9i, 


As : 


ait) 


= 0.07 + 0.92i. 
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final state 


18^/F 


a = 5,/3 = l,5 = l 


5 = 


PP 


8a + 10/3 


50 


50 


nn 


2a + 10/3 


20 


20 




8q + (8 + 25)P 


50 


48 




5a + (5 + 25)P 


32 


30 




2a+ (2 + 2(5)/3 


14 


12 


AA 


(5 + 45)a/3 + (25 + 2(5)/3/3 


24 


16.67 




25a + (8 + 25)13 


20 


8 




25a + (2 + 25)13 


14 


2 




(«-/3)\/3 


6.93 


6.93 



Table 4. The limiting behaviour of 77 —> BB amplitudes. 



Na and Nfj are signature even whereas A^-^ and N5 are sig- 
nature odd. N are isospin 1/2 and A are isospin 3/2, as 
usual. Both exchanges are allowed, although A exchange 
can only take place in pp pp and not in the other sub- 
processes. 

In principle, we should include all four contributions. 
In practice, however, we found phenomenologically that 
the inclusion of just one trajectory, the trajectory, is 
sufficient. 

For the explicit a' in eqn. ([55)1 . as opposed to the a' 
implicit in the trajectory a{t), we adopt 0.9 GeV~^. 

In fig. [21 we show the scaling amplitude, the Regge 
amplitude, and the sum of the two. We first fix the nor- 
malization of the scaling contribution by fitting by the 
eye with the data at W^^ near 4 GeV. We then adjust the 
Regge contribution so that the sum of the two terms fits 
the integrated cross section. 



a(77 — > pp)/nb, |cos&*|<0.6 




0.001 I • • i — • i LI 

2 2.5 3 3.5 4 

W-n /GeY 



Fig. 3. The 77 ^ pp integrated cross section in the region 
I cos 0*1 < 0.6. We show the Belle result against three theoreti- 
cal results: scaling based on the quark-counting rule, the Regge 
amplitude, and the sum of the two. The vertical error-bars only 
includes the statistical uncertainty. 

There is good agreement with the data, except in the 
region just above the threshold. Even this region shows 



improvement compared with our previous calculation in 
ref. [H], where the signature term was neglected. We note 
that by further modifying the signature term by the ar- 
tificial substitution r — z, we were able to obtain the 
fall-off near threshold seen in the real data. This sug- 
gests the possibility that the inclusion of other trajectories 
and/or resonances with appropriate strengths may change 
the threshold behaviour. 

The angular distributions are shown in fig.lH The shift- 
ing of the peak of the central angular distribution from the 
cos 6* « region to the forward region occurs slightly faster 
(about 100 MeV faster) in the theoretical curve than in 
the experimental data. However, the overall trend is in fair 
agreement with that seen in the experiment. 

Having achieved this level of agreement, it becomes de- 
sirable to be able to extend our results to the case of other 
baryons, for instance A and S'^ |24ll25j. However, there is 
no good method for estimating the Regge couphngs [T3]. 
On the other hand, we expect that for all baryons, the 
Regge contributions dominate over the scaling contribu- 
tion, since the Regge contribution is expected to be more 
insensitive to the type of the baryon [15] , whereas the scal- 
ing contribution, from tab. [4l is always smaller than the 
proton pair case. 

As noted in ref. [H], just above the threshold region, 
we expect invariance under u <-!■ c?, as opposed to the 
d ^ s symmetry that follows from the perturbative ap- 
proaches |131I26| . This implies, in the s-channel picture, 
the suppression of isovector a component. Hence S^A pro- 
duction would be suppressed. 

The size of the scaling contribution to each baryon 
pair can be estimated from tab. [4l but this is, as seen in 
the above results, small. The more central, or the higher- 
energy region, is expected to have a more short-distance 
character and so the argument of sec. [2] can be applied. 



3.2 Meson pair production 

We now consider meson pair photoproduction. In eqn. (j54[) , 
we set K = 6. In eqn. ([55| . we set £ = 1. 

The signature r in eqn. (j55p can be ±1 depending on 
the spin of the trajectory. However, unlike in the baryonic 
case discussed above, we have spin-degenerate trajectories 
with T = +1 and r = — 1.t = +1 corresponds to the 
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2.30 <W„< 2.40 GeV 2.40 < W„ < 2.50 GeV 2.50 < W„ < 2.60 GeV 




0.1 0.2 0.3 0.4 0.5 0.6 0.1 0.2 0.3 0.4 0.5 0.6 0.1 0.2 0.3 0.4 0.5 0.6 



2.60 < W„ < 2.70 GeV 2.70 < W„ < 2.80 GeV 2.80 < W„ < 2.90 GeV 




0.1 0.2 0.3 0.4 0.5 0.6 0.1 0.2 0.3 0.4 0.5 0.6 0.1 0.2 0.3 0.4 0.5 0.6 



0.18 
0.16 
0.14 
0.12 
0.1 
0.08 
0.06 
0.04 
0.02 






Belle ' — ' — ' 




Scaling / 




Regge - - - -[/ 




Sum — — 1/ 





0.1 0.2 0.3 0.4 0.5 0.6 



0.03 
0.025 

0.02 
0.015 

0.01 
0.005 




Belle 
Scaling 
Regge 
Sum 



0.1 0.2 0.3 0.4 0.5 0.6 



Fig. 4. The angular distribution of 77 ^ pp. The results of three theoretical model calculations are compared against the 
experimental data from ref. JL^. 



exchange of even-spin mesons: 

A(xr{l- a{t)) (1 + cxp(-i7ra(t))) {a' s)"'-^\ (61) 

The odd-spin mesons have t — — I and: 

cx r(l - a{t)) (1 - exp(-i7ra(t))) (a's)"**). (62) 

Adding together the contributions of degenerate trajec- 
tories, the contribution of the signature term in general 
tends to cancel. In the limiting case of perfect cancella- 
tion, there are two possibilities: 

1. The two amplitude add with the same sign, leading to 
the cancellation of exp{—i7ra{t)). This corresponds to 
the cancellation of the handbag diagrams in fig. [51 

2. The two amplitudes add with the opposite sign, lead- 
ing to the cancellation of the constant-phase term. This 
corresponds to the the cancellation of the cat's-ears di- 
agram in fig. [21 



Let us denote these respectively as 'Regge cat's-ears' and 
'Regge handbag'. It turns out that the plateau structure 
of the K^K^ integrated cross section is only reproduced 
in the 'Regge handbag' case, since a rotating phase is nec- 
essary to yield non-trivial interference with the scaling 
contribution. 

Let us therefore consider the 'Regge handbag' case. 
Here, we can write the combined amplitude as a s — t 
dual amplitude. Using the simple Veneziano amplitude of 
ref. [2lj, we are able to simulate both the resonance and 
the pole regions with the expression: 



A cx 



r{l-a{t))r{l~a{s)) 

ril 



+ (w ^ t), 



(63) 



a{t) — a.{s)) 

from which we can recover the Regge amplitude by the ap- 
plication of the Stirling factorial approximation. Because 
of the resonance-pole duality, the discussion of the pre- 
ceding sec. [21 holds. This has the implication that rela- 
tions between amplitudes such as eqns. (|^ and (|^ are 
satisfied. 
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We fit the scaling amplitude by the eye to the data at 
near 4 GeV. We then study the behaviour of the sum of 
this amplitude and the parametrization of ref. |21j . The 
result is shown in fig. \5\ The scaling curve fits the data 



cr(77^ K+K-) /nh, |cos6i*| < 0.6 



0.1 





Belle ' — ' — ' 




Scaling 


■ H 


Veneziano 




Sum 











1.5 2 2.5 3 3.5 4 

W^^/GeV. 



Fig. 5. 77 K'^ K~ integrated cross section in the region 
I cos 0*1 < 0.6. The experimental results [2] are compared 
against scaling, the Veneziano model, and the sum of the two. 



reasonably above about 2.5 GeV. Below 2.5 GeV, the be- 
haviour of the integrated cross section is still close to the 
scaling curve, but this is accidental since the behaviour of 
the angular distribution is far from that parametrized by 
eqn. ([5^ . On the other hand, the Veneziano amplitude by 
itself provides a semi-quantitative description of the data 
below about 2 GeV although the plateaus just below and 
above 2 GeV are not reproduced. 

The sum of the two amplitudes shows striking resem- 
blance to the real data, except below 1.5 GeV. In par- 
ticular, this reproduces the plateaus. These come from 
the non-trivial interference between the long-distance and 
short-distance contributions to the amplitude, and are, as 
seen in fig. [3 not correlated directly with the shape of the 
resonances. 

This apparent resemblance with the data is, however, 
misleading, since the angular distribution does not cor- 
rectly reproduce the structure of the experimental data. 
This is seen in fig. [SI which shows the distribution at three 
representative energy ranges and in the region between 
2.00 and 2.20 GeV. For most of the energy range, the sum 
of the two amplitudes does not yield a better approxi- 
mation to the angular distribution than either of the two 
individual contributions. 

The plateau structure is weakened but still visible when 
we adopt the Regge limiting expression of the Veneziano 
amplitude, as shown in fig. [T] The situation with respect 
to the angular distribution improves in the Regge limit- 
ing case, especially above 2 GeV. This is shown in fig. [51 
The improvement in the fit compared with the Veneziano 
amplitude is mostly due to the reduced size of the Regge 



cr(77^ if+A'-)/nb, icos6l*j < 0.6 





Regge+scaling 

Belle ' — ' — ' 


- It 













1.5 2 2.5 3 3.5 4 

W^^/GeV. 



Fig. 7. 77 K^K~ integrated cross section in the region 
j cos 0*1 < 0.6. The experimental results [2] are compared 
against the sum of the Regge and scaling amplitudes. 

limiting expression compared to the Veneziano amplitude. 
If the Regge amplitude normalization is modified to more 
closely resemble that of the Veneziano amplitude, the fit 
deteriorates. 

The angular distribution below 2 GeV is not well re- 
produced. 

Since the Veneziano amplitude by itself fits the dis- 
tribution well in the low-energy region, it is tempting to 
introduce a form factor that allows smooth transition be- 
tween the long-distance and scaling amplitudes. We have 
experimented with several such possibilities and found 
that although the fit with the angular distribution im- 
proves, it is difficult to obtain the plateau structure of the 
integrated cross section. 

Let us now turn our attention to the KsKs cross sec- 
tion. We make two assumptions for the K'^K^ = KsKs + 
KlKl cross section with respect to the K^K~ cross sec- 
tion, namely: 

1. the scaling amplitude is scaled by a constant factor of 

1/4; 

2. the Veneziano amplitude remains the same. 

The first implies the complete exclusion of the strange- 
quark contribution, and the second implies the invariance 
of the long-distance amplitude with respect to isospin. 

The results are shown in fig. [9] The large ratio between 
K~^K~ and KsKs cross sections, as well as the increasing 
ratio between the two, is reproduced. However, the ratio 
seems to increase more rapidly in the real data. From our 
above results of the K^K^ cross section, we expect that 
the angular distribution may not be reproduced correctly. 
However, if the long-distance effects are still active in the 
high-energy range, a model-independent statement is that 
the angular distribution would be affected, and will not be 
given by a simple scaling form such as eqn. (|54p . 

We note that for pions, both in the analysis of sec. [51 
and in the Regge/ Veneziano amplitudes, isospin invari- 
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1.40< W„< 1.44GeV 1.80 <W^< 1.84 GeV 2.20 < W„ < 2.24 GeV 




Fig. 6. The angular distribution of 77 — » K~^K~ at three representative energy ranges in between 1.40 and 2.24 GeV (upper 
row), and between 2.00 and 2.20 GeV (lower row). We show the Belle data, the scaling contribution, the Veneziano- model 
contribution and the squared sum of the two. The vertical error-bars on the Belle data are statistical only. 



1.40 <W„< 1.44 GeV 1.80 <W„< 1.84 GeV 2.20 < W„ < 2.24 GeV 




I ' ' ' ' - - T - , - 1 r " " ' - - ' ' 1 

0.1 0.2 0.3 0.4 0.5 0.6 0.1 0.2 0.3 0.4 0.5 0.6 0.1 0.2 0.3 0.4 0.5 0.6 



2.00 < W.„ < 2.04 GeV 2.08 < W.„< 2.12 GeV 2.16 < W„< 2.20 GeV 




I - - - , - - : 1 I ' ' - - " ' 1 ' " ' " - ' ' ' ' 1 

0.1 0.2 0.3 0.4 0.5 0.6 0.1 0.2 0.3 0.4 0.5 0.6 0.1 0.2 0.3 0.4 0.5 0.6 



Fig. 8. The angular distribution of 77 K^K^ , at three representative energy ranges in between 1.40 and 2.24 GeV (upper 
row), and between 2.00 and 2.20 GeV (lower row). We show the Belle data against the addition of the scaling and Regge 
amplitudes. The vertical error-bars on the Belle data are statistical only. 



ance is respected. Hence we expect: 

apj-^ = i (64) 

always, so long as the small difference in the neutral and 
charged pion masses can be neglected. The ratio would be 
violated by the inclusion of 'cat's ears'-type diagrams. 

4 Conclusions 

We studied exclusive hadron pair photo-production pro- 
cesses in low-energy photon-photon collision. 



Motivated by the experimental observation of the un- 
expectedly large suppression of the ratio KsKs/K^K~ , 
we looked into the SU{3) structure of the couplings in- 
volved in these processes, adopting an .s-channci picture. 

We presented the calculations both for the nonet mesons 
and for the octet baryons. We argued that the difference 
between KsKs and K~^K~ cross sections is due to the si- 
multaneous presence of ///' and a in the s-channel. The 
ratio is further enhanced when the /' contribution is sup- 
pressed. We argued that this could be as large as 32. 

We proceeded with a model in which the ratio applies 
predominantly to the part of the amplitude which obeys 
a scaling behaviour. We added to this a long-distance am- 
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a(77 ^ KK) /nb, ] cos 9* \ < 0.6 




0.01 - 



0.001 r 



2.4 2.6 2.8 3 3.2 3.4 3.6 3.8 4 

W^^/GeV. 

Fig. 9. 77 ^ KK integrated cross section in the region 
I cos 0*1 < 0.6 , for K~^K~ and KsKs- We show the region 
2.4 GeV< Wj-, < 4 GeV. We show the experimental results [2] 
and the sum of the scaling and Veneziano model amplitudes. 



plitude whose limiting behaviour is given by Regge theory. 
This latter amplitude tends to be invariant under isospin. 

For the 77 pp process, adopting a Regge amplitude 
for the long-distance dynamics, we obtain a distribution 
that not only fits the integrated cross section well but 
reproduces the behaviour of the angular distribution. 

For 77 —f K~^K~, the long-distance dynamics was 
simulated using the Veneziano model. The summation of 
this and the scaling amplitude results in a curve for the 
integrated cross section that is similar to the experimental 
data. On the other hand, the angular distribution is not 
well reproduced. The Regge limiting amplitude has better 
fit with the data, mainly because of the reduced size of 
the amplitude compared to the Veneziano amplitude. In 
any case, the distribution near the /' peak is reproduced 
better by the pure Veneziano/Regge amplitudes, possibly 
indicating the incompleteness in our parametrization of 
the resonance region. 
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